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ABSTRACT 

Bright emission nebulae, or HII regions, around hot stars are readily seen in Ha light. However, the 
all-pervasive faint Ha emission has only recently been detected and mapped over the whole sky. Mostly 
the Ha emission observed along a line of sight is produced by ionised gas in situ. There are, however, 
cases where all or most of the Ha radiation is due to scattering by electrons or dust particles which are 
illuminated by an Ha emitting source off the line of sight. Here we demonstrate that diffuse, translucent 
and dark dust clouds at high galactic latitudes are in many cases observed to have an excess of diffuse 
Ha surface brightness, i.e. they are brighter than the surrounding sky. We show that the majority of 
this excess surface brightness can be understood as light scattered off the interstellar dust grains. The 
source of incident photons is the general Galactic Ha background radiation impinging on the dust clouds 
from all over the sky. 

Subject headings: ISM: general; clouds; dust, extinction; HII regions; radiative transfer; scattering 
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1. INTRODUCTION 

Although in most instances the Ha emission observed 
along a line of sight is produced by ionised gas in situ 
there are cases where part or all of it is due to scattering. 
In our Milky Way the Ha line radiation from the bright 
ionised regions in the mid-plane of the Galaxy is scattered 
off the widely distributed dust and gives rise to a diffuse 
Ha intensit y rfTura lll979tlWood fc Revnoldsl ll999D. Early 
attempts bv lRevnolds et al.l (|1973| ) to detect the scattered 
Ha light in high latitude dust clouds have not been suc- 
cessful. 

Localised interstellar dust clouds, especially small size 
globules, are frequently found as dar k patches against th e 
bright Ha background in HII regions (|Bok & Reillvill947fl . 
However, the opposite phenomenon of "bright dark neb- 
ulae" should also be observable because of the scattered 
Ha light. This situation is expected to occur when a dust 
cloud is projected towards a faint Ha background at high 
galactic l atitudes. In this letter we show that the all-sky 
Ha map (|Finkbeinerll2003h reveals the presence of numer- 
ous such clouds. Depending on the location of the ionised 
Ha emitting gas along the line of sight, i.e. the proportion 
of gas in front of vs. behind the dust cloud, the cloud may 
actually appear brighter, darker or neutral relative to its 
surrounding Ha surface brightness. 

2. OBSERVATIONAL DATA 



Using the all-sky Ha map (jFinkbeinerl I2003T ) we have 
surveyed a sample of high and inte rmediate galact i c lati- 
tud e clouds listed in th e catalogs of iMagnani et a l. (1985) 
and lReach et al.l (|1998t ) and visible in IRAS 100 //m maps. 
We have found several clouds with an excess Ha sur- 
face brightness but also cases where the cloud appears 
in absorption or even cases with neither "emission" nor 
"absorption" even though the dust column density would 
predict a clear effect. A list of high galactic latitude 
clouds with an Ha excess probably caused by scattering 
includes the following objects: the Lynds clouds L 79, 
L 134/169/183 complex, L 1407, L 1495, L 1532, L 1660, 



L 1780, the MBM objects (jMagnani et al.lll985l ) MBM 16, 
MBM 4 6, MBM 137 MBM 145, and the diffuse infrared 
clouds (|Reach et al.l [l998) DIR 002+31, DIR 046+37, 
DIR 340-43. In Figure 1 we display the Ha and 100 /zm 
images for a selection of four high galactic latitude clouds 
together with the scatter diagrams of Ha intensity vs. 100 
(im intensity. 

The detected excess emission is in the range of 1 to 3 
rayleigh (R) which, when compared with the average all- 
sky Ha surface brightness of ~8 R, suggests that scattering 
can be the source of this phenomenon. The morphological 
appearance of the Ha intensity distribution is very similar 
to the dust distribution as revealed by the IRAS 100 /jm 
images (Fig. 1) and optical extinction maps. The scatter 
plots oi lna vs. I(lOO^im) (Fig. 1 right panel) and Ihu vs. 
A v (Fig. 2 left panel) show for the clouds L134/L183 and 
LI 780, which have low background Ha intensity, a func- 
tional dependence characteristic of scattered light: first a 
linear increase in the optically thin domain, Ay < 1 mag, 
and a gradual saturation for the larger optical thicknesses. 
Other clouds, appearing in "absorption" rather than in 
"emission" (Fig. 1 d; Fig. 2 c-e), can also be readily un- 
derstood as combination of Ha scattering and extinction 
of Ha emission from beyond the cloud's distance. 

The intensity of the Ha radiation along the line of sight 
towards an obscuring and scattering dust cloud at a dis- 
tance d, with the optical depth r at 6563 A , and with a 
geometrical depth through the cloud much smaller than 
the distance d, can be expressed as 

Ii a = Iua (r<d)+ 7 H a (r > d) x e~ T + 1^ (1) 
whereas in the surroundings of the cloud the intensity is 

JBS* = fe« (r < d) + I Ba (r > d) (2) 
and the difference is 

A/Ha = ~ / H a r = S ~ Wr > d) X (1 - O (3) 

Thus, depending on the absolute values of the two terms 
on the right hand side the cloud's Ha surface brightness 
can be higher, lower, or equal to its surroundings. 

The concept of Ha scattering is tested by comparing the 
observed Ha surface brightness with predictions of radia- 
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tive transfer models. We consider a spherical model cloud 
which is illuminated by the Ha Interstellar Radiation Field 
(ISRF). 

3. RADIATIVE TRANSFER MODELING 

The simulations of scattered Ha in tensity were per- 
forme d with Monte C arlo me t hods (|Juvela fc Padoanl 
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pling of the outcoming scattered radiation. The intensity 
of the background sky illuminating t he model cloud was 
obtained from the WHAM survey (jHaffner et alj I2003D 
as function of radial velocity. The missing southern sky 
was filled in by mirroring the data with respect to the 
Galactic longitude. The density distributions of the spher- 
ically symmetric model clouds were calculated according 
to the Bonnor-Ebert solution (|Bonnorlll956l: lEberdri957h 
of hydrostatic equilibrium with temperature 12 K and the 
nondimcnsional radius parameter £ = 7.0. The sizes and 
masses of the clouds were fixed by the visual extinction 
through the cloud centre. Dust properties were based on 
Draine's (2003) 'Milky Way Dust' model. 

During each run, photons of a single frequency were 
simulated and the scattered intensity, including multiple 
scatterings, was registered as an image of the cloud to- 
ward the observer. The full Ha profile was covered at 
steps of 5km/s. The direction of the observer was deter- 
mined by the actual galactic coordinates of the cloud in 
question. This way, although the model clouds themselves 
are spherically symmetric, the anisotropic radiation field 
can produce asymmetrical surface brightness distributions. 
The resulting images of scattered Ha intensity consist of 
21x21 independent pixels. The Monte Carlo noise, i.e., 
the rms-variation of pixel values is about 1%. 

The WHAM survey gives the distribution of the Ha 
emission as seen from the Earth. Although our sources 
are relatively nearby objects, the distribution of the sky 
brightness, as seen from their location, can be different. In 
particular, part of the Ha emission seen toward a cloud at 
a distance of d may, in fact, reside between the cloud and 
us. Therefore, we have also examined scenarios, where a 
given fraction of the Ha emission seen in an area of the 
sky toward and around the cloud is assumed to emanate in 
front of it. Thus, the Ha intensity is divided into the com- 
ponents Iifa(r < d) and I Hair > d) (see Eqs. 1 and 2). 
As seen from the cloud's location the emission component 
I Ha (r < d) is moved to the opposite part of the sky. This 
modification is applied to a region within 45 degrees of the 
direction towards the cloud where the fraction Ih a ( r < d) 
is subtracted from Ifj^" '. The modification is taken into 
account as an additional light source in the opposite part 
of the sky when simulating the scattering of Ha radia- 
tion. The surface brightness component lHa( r < d) is n °t 
attenuated by extinction in the cloud. 

The Ha surface brightness behind the cloud, Isr Q (r > 
d), was originally set equal to the value observed in the 
direction of the cloud in the WHAM-survey. However, the 
resolution of that survey is only one degree, and the surface 
brightness can have large variations within each WHAM 
beam. Of course, the WHAM value also includes the ob- 
ject under consideration whereas in the modeling we need 
a value for the background without the foreground source. 



Therefore, the background sky brightness was adjusted so 
that its level corresponded to the bac kground value ob - 
served in the higher resolution Ha map (Finkb emer [2003h . 

We examine in more detail the Ha surface brightness 
of the clouds L 1780, L 1642, MBM 105, MBM 128/129, 
and MBM 144. The clouds were modeled as Bonnor-Ebert 
spheres such that the range of extinction corresponded to 
the observed values. There is only one free parameter that 
can be adjusted in order to find the best correspondence 
between observed and modeled Ha intensities, namely the 
fraction fcbg, < fcb g < 1, of the background Ha inten- 
sity that resides behind the source, i.e. Jjj a (r > d) = 
fcbg^Ha(surr). The remainder, (1 — fcb g )-Tffa(surr), comes 
from the medium between the cloud and the observer. Fig- 
ure 2 shows the observed and modeled surface brightness 
as function of visual extinction for four clouds. 

In the case of L 1780 the surface brightness can be ex- 
jlained purely as scattered radiation. For the dust model 
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(|Draindl2003[ ) used here the albedo is a — 0.67 and the 
scattering function is moderately forward directed with 
an asymmetry parameter g = 0.5 at the wavelength of the 
Ha line. In L 1642 neither an Ha excess nor absorption 
was observed. This can be explained by means of a suit- 
able fraction of the Ha emission coming from the medium 
behind the cloud. In that case the scattered light from the 
cloud is compensated by the attenuation of the background 
emission, Ihu{t > d). As a result, the surface brightness 
variations remain very small. With fcb g = 0.4 the model 
predicts a total variation of only 1R (see Fig. 2b). In 
MBM 105, because of the higher background level, the in- 
tensity decreases with Ay. The intensity towards the cloud 
centre is determined almost purely by the scattering. In 
MBM 128/129 and MBM 144, with an even higher back- 
ground intensity, the importance of scattering decreases 
and the observations trace mostly the extinguished back- 
ground Ha emission. When a value fcb g = 0.7 is used, i.e., 
when some 30% of the diffuse Ha emission is assumed to 
originate in front of the clouds, the models are in good 
agreement with the observations (see Figs. 2c and 2d). 

Our model clouds are spherically symmetric. However, 
the impinging Ha radiation field is anisotropic. This 
means that the intensity of the scattered light, as observed 
over the cloud's face, is no longer spherically symmet- 
ric. The scattered light was registered as two-dimensional 
maps over the cloud's face. As an example, we show in Fig- 
ure 3 the calculated surface brightness maps for L 1780. 
The asymmetry is quite noticeable at Ay — 4 m . The 
lower frames show predictions for the actual total surface 
brightness where the attenuated background component 
has been added. Although this component is spherically 
symmetric it accentuates the asymmetry of the scattered 
light by reducing the total range of surface brightness val- 
ues. The intensity maximum is clearly displaced towards 
the south, i.e. towards the stronger radiation coming from 
the direction of the Galactic plane. 

4. DISCUSSION 

As shown in Figure 2 a the val ues of a = 0.67 and 
g = 0.5, valid for the dust model (|Draind [2003T ) used in 
this study give a good agreement between the observed and 
calculated scattered light intensities for L 1780. The differ- 
ence can be easily explained, apart from a higher albedo 
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value, by a slightly higher intensity of the Ha radiation 
field at the location of the cloud, a value fcbg < 1, or in- 
homogeneity of the cloud that would in the cloud centre 
cause a higher ratio between scattered intensity and the 
line-of-sight optical depth. For a given cloud direction (I, b) 
a range of (a, g) combinations can reproduce the observed 
scattered light intensity. 

The asymmetry of the Ha intensity distribution is read- 
ily seen in the observations of L 1780 and L 183: the Ha 
intensity maximum is displaced toward the Galactic plane 
relative to the peak of the dust column density maximum 
(see Figs. la,b). This shift is ^8 arc min in the case of 
L 1780. 

No other mechanisms except scattering are needed to 
explain the faint Ha excesses of 1-3 R as observed in high 
latitude translucent dust clouds. The high latitude clouds 
are exposed to the ionising Lyman continuum photons 
which escape from HII regions around O and early B type 
stars near the Galactic plane. The Ha in situ emission 
expected from off-the-plane clouds has been estimated us- 
ing empirical data for the dis tribution of hot stars and the 
intervening dust extinction (Bland-Hawtho rn et al.lfl 998: 
Bla nd-Hawthorn fc Malonevl 1 1999al l2002t iPutman et all 
2003:). For nearby clouds, \z\ < 500 pc, the expected Ha 
surface brightness is < 0.3 R. On the cloud surface facing 
the Galactic plane a thin layer of ionised hydrogen will 
form. Its thickness is As w 5 10~ 2 pc/?i# which for an 
assumed hydrogen density in the cloud's outer layers of 
tih = 10 — 100 cm~ 3 corresponds to 5 10~ 4 - 5 10 -3 pc. 
Morphologically, such a thin H + surface layer would be ob- 
served on the side facing the Galactic plane as uniformly 
distributed surface emission with no correlation with the 
dust column density, and as a thin bright rim around the 
cloud's edge. However, the angular thickness of such a 
bright rim would be of the order 1-10 arc seconds and 



could not be discerned by the present lower resolution data 
for L 1780, the L 134 complex (Figs. la,b), and other sim- 
ilar clouds. 

The conclusion by Idel Burgo fc Cambresvl (|2006f ) that 
the main mechanism for producing the Ha excess in L 1780 
is the ionisation by an enhanced cosmic ray flux is not sup- 
ported by our findings. However, some part of the emission 
could still come from such a process. 

An important aspect of the scattered Ha line radiation 
as compared with the continuum radiation is the velocity 
information it carries both on the source of the imping- 
ing ISRF light as well as on the sca ttering dust particles 
(jBland-Hawthorn fc Malonevl Il999bh . Figure 4 shows an 
example of the observed Ha scattered light line profile. 
L 1780 shows an excess line profile centered at Vlsr ~ 1 
km/s which differs from the mean velocity of the in situ 
Ha emission in the surroundings of the cloud, Vlsr ~ -7 
km/s. This is as expected since (1) the length of the ve- 
locity vector of the cloud itself is probably not much larger 
than its radial velocity component, Vlsr — 3.5 km/s, and 
(2) the impinging ISRF radiation as integrated over the 
whole sky, has a mean velocity of ~ -1 km/s. Model cal- 
culations as function of velocity of the impinging Ha radi- 
ation have been performed for L 1780. We have adopted 
three different scattering functions ranging from isotropic 
(g = 0) to a strongly forward directed one (g — 0.9). While 
all the calculated scattered light line profiles have closely 
the same mean velocity the line width is smaller for the 
forward scattering as compared to the isotropic scattering 
case. In the latter case the cloud "sees" with equal weight- 
ing the whole sky, with impinging Ha velocities ranging 
from Vlsr ~ -50 to +50 km/s. The differences are, how- 
ever, too small to allow a distinction between the different 
values of g (see Fig. 4). 
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Fig. 1. — Fig. 1 Ha and 100 fim images and I Ha vs. I(lOOfim) scatter plots towards some high-latitude molecular clouds. The leftmost 
column shows He* intensity in rayleighs (R) and the central column the 100 fim surface brightness in MJysr - 1 . The intensity scales are 
indicated by the clour bars on the right hand side of the images. The rightm ost column shows pixel-to-pixel rel ation between H c* and 
100 fim intensities. IRAS data have been extracted from the all-sky IRIS map (jMivillc-Dcschcncs & Lagachc 19 93). The pixel size in the 
images is 3.44x3.44 arcminutes, which gives a pixel area that is equal to the native pixel area of the HEALPhc lGorski et ai][20Ml data 
used to make the maps. The native resolution of the He* data is 6 arcminutes. The resolution of the IRAS 100 fim data has been convolved 
into the same resolution. The axes are given in galactic coordinates. Notes on individual clouds: (a) Lynds 134 complex consists of 
the clouds L183, L169 and L 134 (from upper left to lower right) The contours show the extinction at _R-band (A ~ 650 nm), derived via 
starcounts Ouve la et al. 2002): 1, 2, and 4 mag. This region is largely free of foreground/background emission of He*, thus the morphological 
correspondence between scattered He* intensity and H-band extinction is good, (b) Lynds 1780 The cloud is seen in He* excess emission 
against a dim H c* background. The contours show the values of 200 fim optical depth ( Ridderstad et al. 2006): 4 10 — 4 , 8 10 -4 , and 12 10 — 4 . 
c) Ly nds 1642 This cloud is seen neither in excess emission nor in absorption of H c*. The contours show the values of 200 fim optical depth 
Lehtinen efahl 12004*): 410 -4 , 810 -4 , and 12 10~ 4 . The group of bright pixels at I ~ 211.3° , b ~ -35.3° in the He* map is an artifact 
caused by incomplete subtraction of bright stars, (d) MBM128/129 This is an example of clouds seen in absorption against a bright He* 
background. The two minima in the H c* map correspond to maxima in the 100 fim map. 
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Fig. 2. — Fig. 2 Observed (left frames) and modeled (right frames) Ho surface brightness as function of visual extinction. For L 1780 (a) 
and MBM 105 (c) the plotted model values correspond to fc bg = 1. For L 1642(b), MBM 128/129(d), and MBM 144(e) results are shown 
for three values of fcb g . The magenta curves show the average Ha vs. Ay dependence of the best fitting model. The red curves show the 
contribution of scattering, and the blue curves the sum of the attenuated background intensity and, in cases with fc bg < 1, the H-a emission 
that is assumed to reside between the cloud and the observer. For L 1780 we also plot cases a = 0.5 or a = 0.8. The scatter of in the calculated 
intensities is caused mostly by the anisotropic illumination which makes the surface brightness distribution asymmetric (see Figure 3). A 
smaller part of the scatter is caused by the Monte Carlo noise of the radiative transfer calculations. The Av values were estimated using the 
reddening of the background stars that wer e observed in the near-inf rared 2MASS survey. The colour excesses of the stars are converted into 
visual extinction using the NICER method | Lpmbardi & Alvcs (2001). The resolution of the extinction map was set equal to the resolution of 
the Ha map, i.e., 6 arc minutes. 
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Ay=2 m Av=4T 




Fig. 3. — Fig. 3 Modeled Ho surface brightness distribution in L 1780. The maximum extinction through the cloud is assumed to be 
either Av = 2 m (frames a and c) or A\j = 4 m (frames b and d). The upper frames show the distribution of the scattered Ha radiation, and 
the lower frames the total surface brightness where the attenuated background intensity has been added (see Eq. 1). The contours are drawn 
at intervals of 0.5 R. The asymmetry of the surface brightness is caused by the anisotropy of the illuminating ISRF. The maps are oriented 
in Galactic coordinates so that the radiation from the plane of the Galaxy comes from below. The cloud size is specified only by its optical 
depth. 
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Fig. 4. — Fig. 4 Observed and modeled Ho spectrum of L 1780. The figure shows the observed spectrum towards L1780 (uppermost black 
line), the average background around the cloud (black dash-dotted line), and their difference (lower black line). The modeled average signal 
over the source (/total) and the background-subtracted spectra (/total — Ibg) are shown for Draine's (2003) dust model (red lines) and for 
modified dust models with isotropic scattering (g = 0, blue lines) or strong forward scattering (g = 0.9, green lines). The thin coloured lines 
indicate the corresponding spectra for scattered light. The calculations are for a visual extinction of Ay = 2 m through the cloud centre. 



